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A comparison between existing nuclear magnetic resonance measurements and calculations of the 
scalar spin-spin interaction (J-coupling) in deuterated molecular hydrogen (HD) yields stringent 
constraints on anomalous spin-dependent potentials between nucleons at the atomic scale (~ 1 A). 
The dimensionless coupling constant g^gp /4tv associated with exchange of pseudoscalar (axion-like) 
bosons between nucleons is constrained to be less than 5 x 10~^ for boson masses in the range of 
5 keV, representing improvement by a factor of 100 over previous constraints. The dimensionless 
coupling constant fl^PA /47r associated with exchange of an axial-vector boson between nucleons is 
constrained to be g^<7;A /^^i" < 2 x 10~^^ for bosons of mass < 1000 eV, improving constraints at this 
distance scale by a factor of 100 for proton-proton couplings and more than 8 orders of magnitude 
for neutron-proton couplings. 

PACS numbers: 14.80.Va, 11.40.Ha, 21.30.-x 
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Over the past few decades, searches for anomalous 
spin-dependent forces have drawn considerable interest as 
a signature of the axion [T] , a hypothetical pseudoscalar 
Goldstone boson [21 13] arising out of the Peccei-Quinn 
solution to the strong-CP problem [HE]. Axions are also 
appealing as a candidate for dark matter [6]. Other ex- 
otic spin-dependent interactions are predicted by a vari- 
ety of novel theories such as those involving para-photons 
[7] and unparticles [Hj .The possible spin-dependent forces 
that could arise from exchange of scalar/pseudoscalar or 
vector gauge bosons are enumerated in Ref. [9]. An- 
other theoretical framework for considering anomalous 
spin-dependent interactions comes from introducing non- 
zero torsion into general relativity, which causes gravity 
to acquire scalar and vector components that manifest 
as new spin-mass and spin-spin couplings [T0Ul3j . Spon- 
taneous Lorentz violation may also generate exotic spin- 
dependent interactions [14) . 

Recent experiments have significantly improved con- 
straints on long-range (tens of cm) dipolc-dipolc inter- 
actions between neutrons [151 116| and monopolc-dipole 
interactions between electrons and nucleons [l7]. Earlier 
work constrained monopole-dipole interactions between 
nuclei [18 and anomalous dipole-dipole interactions be- 
tween electrons [THl [20] and electrons and nuclei [H]. 
Constraints on new pseudoscalar interactions can also 
be obtained from spin-independent tests of the inverse 
square law [22] • Other laboratory experiments constrain 
anomalous monopole-dipole couplings on length scales in 
the range of /im to mm [23H26J . Spin relaxation stud- 
ies of hyperpolarized '^He gas have been used to limit 
anomalous dipole-dipole interactions between neutrons 
at length scales of about 100 nm ^7\. Because of the 
possibility that new force-mediating bosons may be mas- 
sive and have limited range, it is of considerable interest 



to find experimental techniques to search for anomalous 
spin-dependent interactions at even shorter distances. 

Here we discuss constraints on the existence of anoma- 
lous dipole-dipole forces on angstrom length scales. Our 
constraints are obtained by comparison of nuclear mag- 
netic resonance (NMR) measurements and theoretical 
calculations of J-coupling in deuterated molecular hydro- 
gen (HD). Such couplings have the form JI- S (here I and 
S are nuclear spin operators) and arise due to a second- 
order hyperfine interaction. To our knowledge, the only 
constraints on anomalous nuclear dipole-dipole couplings 
at these length scales come from Ramsey's molecular- 
beam measurements of H2 [28] and studies of spin- 
exchange collisions between "^He and Na [29]. Atomic- 
scale constraints on anomalous, dipole-dipole couplings 
between electrons and nuclei, mediated by exchange of 
axial-vector bosons, have recently been obtained from 
the hyperfine structure of hydrogen- like atoms [5nH32[ . 
The constraints derived in this letter represent a factor 
of 10^ to 10^ improvement over previously derived labo- 
ratory limits on anomalous spin-dependent forces at the 
A range, corresponding to a particle with keV-scale mass. 

There are several additional noteworthy features of the 
analysis presented here. Constraining monopole-dipole 
interactions is appealing because such couplings violate 
invariance under both time reversal (T) and spatial in- 
version (P) , and hence one expects negligible background 
from standard model physics. Dipole-dipole couplings, 
on the other hand, are even under both T and P, and 
arise from standard model physics. In this sense, dipole- 
dipole couplings may appear less attractive in searches for 
exotic physics because one must carefully calculate the ef- 
fects of standard model physics. Despite this difficulty, 
our analysis results in constraints that are within two or- 
ders of magnitude of hadronic axion models. This is sig- 



nificantly closer than limits on monopole-dipole couplings 
come to constraining QCD axions [23l - f26] . Furthermore, 
axion mediated dipole-dipole coupling scales as !//„, {fa 
is the Peccei-Quinn symmetry breaking parameter) inde- 
pendent of the QCD 9 parameter. In contrast monopole- 
dipole and monopole-monopole coupling depend on 9 as 
9/fj^ and 9^/f^, respectively [T]. Therefore, interpreting 
limits on monopole-dipole or monopole-monopole cou- 
plings as a constraint on fa comes with large uncertainty 
related to the size of the 9 parameter. In many discus- 
sions, 9 is taken to be order of 10"^" based on limits from 
the neutron electric dipole moment experiments, however 
it may be much smaller, relaxing the constraints one can 
place on fa- 

We also note that stringent constraints on axion cou- 
plings can be obtained from astrophysical observations. 



Electron-axion or photon-axion couplings can be con- 
strained by cooling rates of low-mass stars [33l |34] , re- 
quiring the axion mass be less than 10~^ eV. More rele- 
vant in the present context are the constraints obtained 
from supernova SN1987A [3S] or the metallicity of red 
giant stars |36| . which limit axion-nucleon interactions. 
The 14.4 keV emission line of ^^Fe also constrains axion- 
nucleon couplings |37]. However, these astrophysical con- 
straints are somewhat particular to the axion, and do not 
apply to other pseudoscalar bosons or axial-vector inter- 
actions PISH]. 

Moody and Wilczek discuss the potentials arising from 
the exchange of pseudoscalar (P) axion-like particles [1] 
in the non-relativistic limit. The dipole-dipole potential 
has the form 
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Here, we work in units where h = c = 1, ?7i is the mass 
of the axion-like particle, gpgp/{4:TT) is the dimensionless 
pseudoscalar coupling constant between the particles, Mi 
and M2 are their respective masses, and r is the distance 
(in units of inverse energy) separating the two particles. 
The subscript on the left-hand side is chosen to match the 
notation of Ref. [9] . We also note that a S function con- 
tribution to V3(r) is neglected here because of Coulomb 
repulsion of the two nuclei. 

The measurements from which we extract our con- 
straints occur in gas phase, in which the internuclear 
vector f suffers random reorientation due to collisions, 
leading to averaging of Eq. (II]) . The first term in braces 
has the same angular dependence as the usual magnetic 
dipolar interaction and averages to zero. After averag- 
ing, the second term is proportional to cti • 0-2 , yielding an 
effective anomalous J-coupling AJ3I • S, where I and S 
are the respective spins of the proton and deuteron, with 
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Here, gp = gp + gp, we have made the approximation 
that the neutron and proton masses are equal, Mn = 
Mp, and we assume that the proton and neutron of the 
deuteron each contribute roughly equally to the spin, S, 
of the deuteron, (dp) = ((T„) = S. 

Two measurements of J in HD can be found in the 
literature: Ref. [Sg reports J = 42.94 ± 0.04 Hz and 
Ref. 40J reports J = 43.11 ± 0.02 Hz at 40 K. The 
mean and standard deviation of these measurements is 
43.025 and 0.12 Hz, respectively. Using density func- 



tional theory, Vahtras et al. |41j calculate a spin-spin 
coupling constant of J = 43.15 Hz at 40 K, in agreement 
with measurements at the level of 0.12 Hz. No estimate 
for the uncertainty of the calculation was reported. We 
take this to constrain AJ3 < 0.24 Hz (9.8 x 10" ^^ eV) 
at the 2-sigma level. The distance between proton and 
deuteron is about r — 1.4 Bohr radii |41j . or in units 
where n = c=l,r = 0.00038 eV'^ 

From the experimental measurements, theoretical cal- 
culations, and Eq. (pi), we can constrain the dimension- 
less coupling parameter gpgp/ATT as a function of the 
pseudoscalar mass m, as indicated by the light shaded 
region bounded by the dashed line in Fig. [T] The 
strongest constraint occurs for bosons with mass m = 
2/r = 5300 eV for which gj^g^p/Air < 5 x 10"''. For com- 
parison, the darker shaded region bounded by the dashed 
line shows the limits obtained from Ramsey's molecu- 
lar beam measurements i28j of HH dipole-dipole interac- 
tions, where non-magnetic contributions are constrained 
to be less than Knm < 70 Hz, (3 x 10"^^ eV), at the 
2 — (T level. The form of the interaction Ramsey con- 
sidered differs from Eq. (IT]) in its angular dependence, 
however it can be approximately interpreted as 
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This limit (2 — tr level) is indicated by the dark shaded 
region bounded by the solid line in Fig. [l] 

We note that our constraints are within two orders of 
magnitude of standard axion models. The axion mass 
is determined by the Peccei-Quinn symmetry breaking 
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FIG. 1: Constraint (at the 2-a level) on the dimension- 
less coupling constant gpg^/A-jr associated with the dipole- 
dipole potential, V3(r), as a function of the exchange boson 
mass, obtained from comparison of measured and calculated 
J-coupling parameters for HD (dashed line, light gray fill), 
where gp ~ gp + G^p. Also shown is the existing constraint 
on p — p pseudoscalar couplings from Ref. |28| (solid line, 
darker fill). The straight fine represents axion couphng in the 
KSVZ model, given by Eq. (|5f . 




FIG. 2: Constraints (at the 2-a level) on the dimensionless 
coupling constants 5a<?a/(4'7I") as a function of the mass of an 
axial vector bosons (A) between nucleons [V2(r), see Eq. (|8|]. 
Curve (1) shows the constraints obtained from the J-coupling 
interaction in HD, where 1 = p and g\ ~ dA '^ 9^- Curves 
(2) and (3) show the respective constraints on e — n and e — p 
axial vector couplings from Ref. [31]. Curve (4) shows the 
constraints obtained from Ramsey's measurements of proton- 
proton dipole-dipole coupling in molecular H2 . 
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where /^ — 92 MeV is the pion decay constant, m^ = 
134 MeV is the pion mass, and z — m^/raj: = 0.56 is the 
ratio of up and down quark masses. In the Kim-Shifman- 
Vainshtain-Zakharov (KSVZ) model (see, e.g. Ref. [H] 
for a brief review) , axions couple only to hadrons but not 
quarks or leptons at tree level. Its coupling to a nucleon 
of mass Mm is given by gp = CNMjy/fa, where Cat is 
a model-dependent numerical factor. Using expressions 
for Cn found in Ref. [32] for KSVZ axions, we find 
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This relation is given by the straight line in Fig. [T] 
Our analysis, obtained from purely low energy laboratory 
measurements, is thus within two orders of magnitude of 
constraining KSVZ axions. 

It is worth discussing axion constraints obtained from 
high energy experiments in the context of our limits. To 
our knowledge, the best limits come from the decay of 
the K~^ particle. Reference [H] limits the branching ra- 
tio BR{K+ ~> Tr+ + a) < 4.5 x 10"". The relation 
between the branching ratio and the Peccei-Quinn sym- 
metry breaking scale can be found in Ref. [H]: 
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Using values of Cn found in Ref. [12], and for z = 
0.56, we find fa > 11600 GeV (corresponding to nia < 
507 eV), and interpreted in the context of Fig. fl] 
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While these high energy limits are significantly stronger 
than those obtained from HD, the latter are specifically 
sensitive to the spin-dependence of the interaction. 

In addition to pseudoscalar axion-like particles, the 
presently considered set of measurements and calcula- 
tions can be used to constrain spin-dependent forces due 
to exchange of other particles. In the case of spin-1 axial 
vector (A) bosons, there arises a Yukawa potential [9] 
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Again assuming that the nucleons in deuterium con- 
tribute equally to its spin, we find that in the limit of 
massless bosons (g^+5^)5^/47r < Jr/2 = 1.8 x 10~^^ at 
the 2 — cr level. Constraints as a function of boson mass 
are shown in Fig. p^ (curve 1), along with constraints ob- 
tained from Ramsey's measurements in H2 (curve 4) . We 
also show the spin-dependent constraints between elec- 
tron and neutron (curve 2) or proton (curve 3), obtained 
from measurements and calculations of the 2s hyperfine 
state in atomic hydrogen and deuterium [3T]. The lim- 
its obtained from spin-exchange cross sections in ^He-Ne 
collisions [29] are about eight orders of magnitude worse 
than obtained in the case of HD J-coupling. 

Finally, we point out that J-coupling measurements in 
an isotropic liquid or gas, where molecules rapidly rotate. 



are most sensitive to exchange of axion-like particles with 
masses in the range of 2/r [see Eq. ^] . If the molecule 
is placed in an anisotropic environment such as a liquid 
crystal, tumbling is inhibited and the molecule becomes 
partially oriented, thereby rendering the term in curly 
brackets in Eq. (fTl) observable. This term is largest for 
pseudoscalar exchange boson mass m = 0. Therefore, 
precise calculations and measurements of dipole-dipole 
couplings in liquid crystal environments may yield im- 
proved constraints on light pseudoscalar bosons. Accu- 
rate measurements and calculations of dipole-dipole cou- 
plings in partially oriented benzene exist |45j . however 
the limits one may extract from them have potentially 
larger systematic uncertainty and will be discussed else- 
where. 

In conclusion, we have used measurements and theo- 
retical calculations of J-coupling in deuterated molecular 
hydrogen to constrain spin-dependent forces due to ex- 
change of exotic pseudoscalar and axial- vector particles. 
This analysis improves by over two orders of magnitude 
the constraints obtained by Ramsey's comparison of ex- 
periment and theory on the dipole-dipole coupling of pro- 
tons in molecular hydrogen. In particular the constraints 
on the exchange of axion-like pseudoscalar particles, is 
improved by a factor of 10^ in the mass range of 10^ — 10* 
eV. We also discussed the limits these measurements can 
place on the exchange of photon-like lightweight axial- 
vector bosons, improving constraints by eight orders 
of magnitude over limits placed on the proton-neutron 
couplings, compared to an analysis of spin-exchange in 
sodium and neon. New experimental techniques to mea- 
sure scalar couplings with extremely high accuracy based 
on ultra-low field NMR [33 could lead to improvements 
in experimental precision. Combined with improved den- 
sity functional theory calculations, comparison of scalar 
or dipole-dipole couplings may be used to further con- 
strain spin-dependent forces. An improvement in the 
pseudoscalar coupling limit by two orders of magnitude 
would bring these constraints into the range of standard 
axion models, providing a new set of purely laboratory 
limits on the QCD axion in the keV range. In order to 
substantially improve these limits it is desirable to find 
spin systems with a naturally small J coupling to reduce 
the requirements for precision calculations of standard 
model interactions. In this regard it is interesting to con- 
sider chemically unbound spin systems. In this case, a 
finite J coupling can still arise due to second-order hy- 
perfine interaction in van der Waals molecules [47]. For 
example, in a mixture of liquid ^^^Xe and pentane, a J 
coupling of 2.7 Hz was measured in agreement with pre- 
dictions [31]. Lighter atoms would have even smaller J 
coupling ^49j. A liquid H2 — '^ He mixture [50] may be 
interesting to consider in this context. Hyperpolarizing 
■^He in the mixture would allow one to measure very small 
frequency shift of H NMR using the techniques described 
in [31 [51]. 
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